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We study theoretically the Raman scattering spectra in the one-dimensional (ID) quantum spin- 
I antiferromagnets. The analysis reveals that their low-energy dynamics is exquisitely sensitive 
to various perturbations to the Heisenberg chain with nearest-neighbor exchange interactions, such 
as magnetic anisotropy, longer-range exchange interactions, and bond dimerization. These weak 
interactions are mainly responsible for the Raman scattering and give rise to different types of 
spectra as functions of frequency, temperature, and external field. In contrast to the Raman spectra 
in higher dimensions in which the two-magnon process is dominant, those in ID antiferromagnets 
provide much richer information on these perturbations. 

PACS numbers: 78.30.-j, 75.10.Pq, 78.67.-n, 75.40.Gb 
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Introduction: Quantum antiferromagnets have long at- 
tracted much attention as a laboratory to study quan- 
tum many-body effects. Experimentally, several tech- 
niques are available to investigate them; measurements 
of magnetic susceptibility, specific heat, and spectra of 
neutron scattering, NMR, and ESR. Recently, the opti- 
cal spectra [ll-Q have also turned out to be a powerful 
tool to study quantum spin dynamics. An example is 
electromagnon spectroscopy of multiferroics, where the 
one-magnon process is activated by an electric field in 
infrared absorption due to the magnetostriction mecha- 
nism [J]. Raman scattering on the other hand, 
is usually considered to detect two-magnon excitations in 
antiferromagnetic (AF) ordered phases, which has been 
utilized to estimate the strength of the exchange interac- 
tion. 

In one-dimensional (ID) systems where quantum fluc- 
tuations are much enhanced, such a simple magnon pic- 
ture fails miserably. A canonical model for ID quantum 
antiferromagnets is the spin-^ Heisenberg Hamiltonian, 
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jY.s,-s, 
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where J (> 0) is the exchange interaction between neigh- 
boring spins. The low-energy physics of this system is de- 
scribed by a Tomonaga-Luttinger (TL) liquid with gap- 
less spinon excitations 0] instead of magnons. In real sys- 
tems, however, additional small perturbations V always 
exist, e.g., spin-orbit interaction, magnetic anisotropy, 
disorder, longer-range exchange interactions, and also 
spin-lattice coupling leading to the bond dimerization 
(spin-Peierls instability). Despite the smallness of these 
interactions, they are crucial for the quantum dynamics 
of the system, and are the subject of intensive studies. 

Unfortunately, experimental signatures of these small 
perturbations V are often difficult to study because con- 
ventional experiments probe quantities that are domi- 



nated by the Heisenberg term ([T]) of the Hamiltonian. 
Therefore, it is highly desirable to have experimental 
probes that reveal the physical processes associated with 
the small perturbations V in ID antiferromagnetic sys- 
tems. In this Letter, we show that Raman scattering 
from ID spin-i antiferromagnets provides such an ex- 
perimental probe. 

It has been considered so far that the Raman scat- 
tering in ID magnets is not so useful compared to other 
conventional methods although some of the experimental 
and theoretical works exist This is because the 

Hamiltonian T-Lq and the corresponding Raman operator 
TIq (X Ho [see Eq. ([3])] commute with each other, and 
hence no Raman scattering occurs without additional in- 
teractions. Furthermore, these perturbations V, which 
determine the Raman scattering spectra (RSS), remain 
rather uncertain in most cases. However, this does not 
necessarily mean that the RSS is useless in these systems. 
In fact, once theoretical predictions on the RSS for each 
interaction V are available, RSS can provide useful infor- 
mation on V as we will see later. 

The results of our analysis based on field-theory and 
nonperturbative methods are summarized in TablesUand 
Fig.[2]for gapless cases, and Table HIl and Fig.[3]for gapped 
cases, respectively. Comparing these predictions with the 
observed temperature, frequency, and magnetic-field de- 
pendence of the RSS, one can obtain detailed information 
on V. The results will be explained below. 

Definition of RSS: Let us start from the definition of 
the RSS and the Raman operator. The RSS is propor- 
tional to the dynamical structure factor of the Raman 
operator TZ, namely. 
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dt 
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where lo — uji ~ ujs and ^^(5) is the energy of incident 
(scattered) photon. In Mott-insulating systems, the Ra- 
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man operator [H, [3, S @ usually has the form 

7^ = ^ (e. • ri2)(e, • ri2)A{ri2)Sr, ■ Sr„ (3) 

ri ,r2 

where 6^(3) is the polarization direction of the incident 
(scattered) photon and ri2 — ri — r^- Therefore, the 
RSS strongly depends on the direction of applied and 
observed electromagnetic waves and the crystal structure 
of magnets. The factor A{ri2) is difficult to accurately 
determine, but the ratio between the factors on different 
bonds is known to be of the same order as that between 
the exchange couplings on those bonds. From Eqs. ([2]) 
and (O, one can easily find that the intensity /(w) is 
unchanged when TZ is replaced with the modified Raman 
operator 

n' = n-cn, (4) 

where C is arbitrary real constant and H is the Hamil- 
tonian of the target magnet. We can therefore adopt TZ' 
to make the calculation of I{uj) easier. 

Analysis: The low-energy physics of the Heisenberg 
chain ([l} with/ without an easy-plane anisotropy Vi — 
— J^YlijSjSj+i s-iid ^ Zeeman term is well described 
by the TL-liquid theory The low-energy effective 
Hamiltonian is identical to the free boson theory 

Hf ^ j dx^-[K-\d,cj,f + K{d,ef] , (5) 

where x = jao (oq is lattice spacing), v is the spinon 
velocity, K is the TL-liquid parameter, and (</>, 9) is the 
canonical pair of bosonic fields. The parameter K = 1 
at the SU(2) symmetric point, while an anisotropy Vi 
or an external field H usually increases the value of 
K, i.e., K > 1. Note that Vi is the perturbation 
in the sense that it violates the commuting property 
of the Hamiltonian Hq + Vi and the Raman operator 
TZ. Spin and dimer operators can also be bosonized 
as Sf « J^'ix) + i-iyWix) and {-lySj ■ Sj+i « 

dsin(v^(/') H , where 7V+ = boe'"^" H , J+ = 

hie"^^ cos(V2^0) + • • • , 7V^ = ai cos(V2^<?!') -f • • • , and 
J'^ = ao9a;(/)/\/27r. The anisotropy and field dependence 
of parameters ai, 6o,i, d, K and v is accurately evalu- 
ated [lB-E3]- The bosonization approach therefore en- 
ables us to estimate the effects of several perturbations 
on the RSS /(w) with reasonable accuracy in the low- 
energy region, i.e., T, \uj\ ^ J. 

Gapless Cases: Let us study four realistic perturba- 
tions V that do not violate the TL-liquid phase; XXZ 
anisotropy Vi, longer-range exchange couplings V2 = 
Y.n>2Y.j JnSj-Sj+n {\Jn\ < J), a boud tilting in Fig. [H 
and a random bond alternation V4 = J2j ' 
Sj+i with Uj being the randomly distributed lattice dis- 
tortion {\uj \ <^ 1). The results are summarized in Table |T] 
and Fig. H 



Generally Vi,2 are always present in real compounds. 
The main part of the bosonized V2 is 

I dx + h.c.) + c^J^Jl (6) 

where is the right (left) moving part of J'", and 

constants Cxy = Cz depend on J and J„. Similarly, we 
obtain Vi w -JA / ^{J'- - Af^Af^). From Eq. Q, 
we can make TZ' proportional to Vi^2 if ^i,s are set parallel 
to rj — Tj+i. Applying the standard technique based on 
the bosonization and conformal field theory f?\, we can 
calculate the Raman intensity of the Heisenberg chain Tio 
with Vi or V2 for arbitrary frequency uj and temperature 
T — 1/(3. The result for the case of V2 is 

/(t^) (X 2clyF{L,,p, K) + clF{i^,p, l + e)U^o (7) 

where F{x,y,z) = ^(l-er^y)-^ sm{2nz)lm[B{~i^ + 
z,l ~ 22)2](2£M)4z-2 j^^^j B{x,y) is the Beta func- 
tion. We also have /(w) oc ^n^ajFiuj, (3, K) + K^i^ + 
^al)^F{uj, (3,1 + e) for the case of Vi. The intensities of 
these two cases are presented in Fig. [3] (a) and (b), and 
they show that /(w) is a monotonically increasing func- 
tion of UJ, and it remains finite in the limit w/J — > at 
finite temperatures. These properties are at least qualita- 
tively consistent with the experimental result of the para- 
magnetic phase of CuGeOa in which J ~ 150K and 
J2 ~ 30K. We note that Vi and V2 give the similar be- 
havior in the RSS, but they can be distinguished by other 
physical quantities such as susceptibilities. When we ap- 
ply a magnetic field H and a magnetization M — {Sj) ap- 
pears, B{-i^ + K,l- 2KY in the first term F{uj,[3, K) 
of /(w) is changed into s(_,Ei£±il^£!V£o) + 1 - 

2jr)g(- /^"-^;f + j^, 1-2K) for both cases of Vi,2. 
As a result, the RSS weight of this term becomes nearly 
zero in the low-frequency region w < wi = AirMv/aQ at 
low temperatures T <^ J. For instance, in the case of V2, 
only the term survives in Eq. ([7]). 

We next consider ID magnets with a tilting bond as 
in Fig. [1] In fact, tilting structures with a small an- 
gle 6*0 exist in several cuprate magnets such as Cu ben- 
zoate ilSi], KCuGaFe [II], and [PM]Cu(N03)2(H20)2 
(PM=pyrimidine) [20[. In this system, the Hamiltonian 
is the same as Eq. ([T]) and hence a TL-liquid state sur- 
vives. However, if we fix e^^g as in Fig. [1] the Raman 
operator becomes different from that of the case with- 
out a tilting bond. Tuning the value of C in Eq. we 
obtain 

TZ' (X sin(26lo) sin(6l, + g^) ^(-l)JS'j • Sj+i. (8) 

j 

This is nothing but a dimerization operator and does not 
commute with Tio- Using this operator, we obtain /(w) cx 
sin2(26io)sin^(6'i + es)d^F{uj, j3,K/2) that is depicted in 
Fig. [21(c) The RSS rapidly increases around w = T 
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TABLE I: Properties of RSS in gapless cases of ID Heisenberg magnet Ho with perturbation V. Constants Ci_i2 depends 
on the TL-liquid parameter K, the spinon velocity v, the magnetization M = (Sj), the lattice spacing ao, etc. The value of K 
is unity at the SU(2)-symmetric case, while an easy-plane XXZ anisotropy or a magnetic field H increases it, i.e., K > 1. 



perturbation V 


bosonized form of V 


scaling dimension 


RSS liu)) 




mam ertcct of held H 


XXZ anisotropy 




2K {N'' term) 

2 ( J"^ and M'" terms) 


cicj*"-^ + C2aj^ (/3a; 
C3r^^-2 + C4r2 (/3c. 


»1) 
«1) 


ci,3 terms disappear 
in cj < oJi = AnMv/aa 


longer-range coupling 

V2 = E, Jr^SJ ■ 5j+„(n > 2) 


C:.,(^^.7^-+h.c.)/2 


2A' (ci:j, term) 
2 (cz term) 




»1) 
«1) 


C5,7 terms disappear 
in cj < oJi 


tilting bond in Fig. [T] 


sin(26lo)sin(6', + 6),) 
xdsin(\/27r<^) 


7^/2 


Cgc"-^ (/3CJ> 1) 




C9,io terms disappear 
m u) < ll)2 = 27rMw/ao 


random dimerization 


Jitjdsin(\/27r<^) 


if/2 


Clio;-"-' (/3a; > 1) 
ci2T^-i (/3cj < 1) 




value of K increases 



^ j j+1 j+2 








e. 



FIG. 1: (color online) ID antiferromagnet with tilting angle 60 
and polarization directions of incident and scattering photon 
Bi.s with angles 6i and 9s. 



at low temperatures, and the form is quite different from 
the case of Vi,2- Physically the origin of this spectrum 
is two-spinon states. We emphasize that the strength of 
can be controlled by tuning angles Oi s- Similarly to 
the case of V1.2, a magnetic field H makes the weight of 
I{uj) absent in the region a; < 0^2 = 2t:Mv / a^. 

A random dimerization is expected to be present in the 
higher-temperature paramagnetic phase of spin-Peierls 
compounds such as CuGeOa [§| and TiOCr [lH. In this 
case, the Raman operator TV is proportional to V4/J ~ 
/ ^ttjdsin(\/27r(/)). Under the assumption that Uj is a 
sufficiently small perturbation from T-Lq and {ujUkjR = 
u'^Sjk {{■■■)r stands for the average over the random- 
ness), the Raman intensity (/(cli))^ is reduced to a local 



correlator cx u^(P J 



It is calculated as 



dt_^iu}t / 
27r ' 



sm{V2n(t>{t)) sin(V2^(/!)(0))). 



v?<Pao / 27rao 
47ri' \ I3v 



K-l 



0^ I K I3uj K I3uj\ , , 



which is shown in Fig. [5] (d). The cj dependence of 
{I{uj))ji is negligible in w > T. Such a spectrum with 
a small slope is observed in the paramagnetic phase of 
CuGeOa (see the region uj < 50cm~^ in Fig. 1 of Ref. @), 
and therefore the spectrum might contain the contribu- 
tion from V4. An applied field H does not affect the 
form of {I{u!))r much, but it slightly varies parameters 
(-ftT, u, ai,6o,i,(i). 

Gapped Cases: Let us now discuss another kind of typ- 
ical perturbations V that break the TL liquid in Hq and 
open a finite excitation gap: a static bond alternation 
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FIG. 2: (color online) RSS /(cj) for ID Heisenberg magnet JT]) 
with additional perturbation Vs. Panels (a), (b), (c) and (d) 
correspond to the cases of Vi, V2, the bond tilting in Fig. [T] 
and V4, respectively. All the continuous spectra come from 
multiple spinon states in the TL liquid phase. 



and 



(dimerization) term V5 = J2j ^i^^Y 
form and staggered Zeeman terms Ve = — J2j HSj , 
{—lyhSj induced by an applied field H. The results are 
summarized in Table HIl and Fig. [S] 

In the case of V5, the effective Hamiltonian becomes 
an exactly solvable sine-Gordon (SG) model. 



njci 



off 



dx 
ao 



-udsm{\/2nq 



(10) 



There are three kinds of massive excitations: soliton (5), 
antisoliton (S), and some breathers (i3„) that are the 
soliton-antisoliton bound states. The mass of the soliton 



Es is equal to that of antisoliton, and it is given by [17 1 



El 
J 



2 



K 



Jqq nd r( 



4 ' 



2 r(x/4) 



(11) 



where T{x) is the Gamma function. The 71-th breather's 
mass En is related to Eg via En = 2Es sin[n7r/(8//'ir — 2)] 
with n = 1, • • • , [4/if — 1]. The SU(2)-symmetric dimer- 
ized chain with K = 1 has only two breathers Si, 2- 
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TABLE II: Properties of RSS /(i^) in gapped cases of ID Heisenberg magnet "Ho with perturbation V at T = 0. In the case 
of Vs, the soliton mass, and first and second breather masses are respectively evaluated as Eg ~ 1.5u^^'^J, Ei = Eg and 
E2 = VSEs, while Es « 2.\{hlJfl'^J in the case of Ve with a small field H <^ J . Note that in the case of V5, E2 becomes 
larger than y/^Es [H if the marginal operator, neglected in the SG model H5, is taken into account. On the other hand, a 
small next-nearest-neighbor AF coupling J2 weakens the effect of the marginal term 



perturbation V 


bosonized form of V 


Raman active mode (its main origin) 


peak positions for each mode 


static dimerization 


Judsin(\/27r<^) 


second breather (sin(\/27r0) term) 
S-S continuum (sin(\/27r<^) term) 


Lj — E2 (stable against H) 
uj > 2Es 


uniform and staggered 
Zeeman terms 


-hbo cos(^2^6l) 


soliton, antisoliton (tilting bond) 
odd-th breathers {dx(l> term) 
even-th breathers {cos{\/2tt9) term) 


CO = {E^ + {2TTMv/a„fy^'^ 

iiJ ~ E2n + 1 
1^ ~ E2n 



Three particles S, S and Bi corresponds to massive spin- 
1 triplet excitations with = +1, —1, and 0, respec- 
tively, while B2 is regarded as a singlet excitation with 
S = 0. The soliton mass is evaluated as Es ~ 3.5((iu)^/'^ J 
with d « 0.3 [13] at the SU(2) point. From Eq. dS]), the 
RSS is proportional to the dynamical structure factor of 
sm{\/2n(j)). To accurately evaluate such dynamical corre- 
lators of the SG model at the low-energy region, we utilize 
the form-factor approach [H, [i^] which is reliable when 
T/J ^ 1. From this approach, the lowest- frequency con- 
tribution of I{uj) is shown to be a 5-functional peak of 
the singlet breather B2 at uj — E2, and the second low- 
est one is given by the continuum spectrum of soliton- 
antisoliton scattering states with oj > 2Es. The weight 
of each contribution can also be exactly calculated by the 
form-factor method. In particular, the weight of the sin- 
glet breather is proportional to {Esao/v)^ , and is much 
larger than that of the continuum. The B2 peak and 
its weight in I{uj) are shown in Fig. [3] (a). The distor- 
tion (u) dependence of this peak can be compared to 
Raman scattering experiments for spin-Peierls magnets, 
CuGeOg TiOCr etc. The B2 peak of /(w) is 

stable against an applied field H ii H is smaller than the 
critical field He — Eg. 

A staggered magnetic field h emerges as we apply a 
uniform field H to magnets with a staggered gyroniag- 
netic tenso r |24l| . Typical examples are Cu benzoate [l8j , 
KCuGaFe n% and [PM]Cu(N03)2(H20)2 [20^, in which 
a tilting structure is also present (as wc discussed). In 
these compounds, h « CgH {\cs\ <^ 1) is realized. The 
bosonized form of Vg is given by 



(12) 



Haodx4'/\'2Tr — hbo cos(v 27r0). 



The term dx(f>, inducing a finite M , can be absorbed into 
the free boson part Hq^ , and then the effective Hamil- 
tonian is also a SG model 2Jj. Therefore, we can again 
apply the form- factor method to calculate I {to). The 
soliton mass is given by Eq. ([TT|) with replacing {K, d) 
by (l/Kjhbo), and the breather masses are given by 
En = 2Es sin[n7r/(8i^ - 2)] with n = 1, • • • , [4/^ - 1]. 
Since the value of K increases with increasing iJ, the 
number of breathers [AK — 1] is also increased with H 




ro/J 

Cs=0.1 
(h=QH) 



Bi 
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FIG. 3: (color online) (a) u dependence of the singlet-breather 
{B2) peak in /(w) of the dimerized magnet Ho -I- V5 at T = 0, 
where 5 functions are broadened. The peak position u) = E2 
and its weight are both proportional to v?^^ . This peak is sta- 
ble against an applied field H. The weight of the continuum 
spectrum with uj > 2Es is much smaller than that of the B2 
peak, and the former is omitted here, (b) H dependence of 
peak positions of each particle in I(io) of the magnet Ho + Ve 
at T = 0. Some level crossings between the soliton <S and 
breather B„ peaks occur with increasing H. The continuum 
spectra are omitted. 



in the present case. From the form- factor method [23|, 
dx4', cos(-\/27r0), and the tilting-bond term sin(-\/27r0) in 
the Raman operator are respectively shown to provide S- 
functional peaks of odd-th breathers at a; = ii'odd, even- 
th breathers at w = -Eovon, and soliton (antisoliton) with 
wavenumber k = 2ti M/a^ at w = [E^ + k'^v^y^^ in the 
spectrum I{uj). Namely, in contrast to the case of V5, 
all of the elementary particles of the SG model can be 
observed. The H dependence of several peak positions 
are plotted in Fig. [3] (b). Remarkably, level crossings 
between soliton and breather peaks occur ^19]. In addi- 
tion to these peaks, there exist continuum spectra with a 
smaller weight, although it is difficult to accurately eval- 
uate them. 

In conclusion, we have shown that various weak per- 
turbations V to the spin-i AF Heisenberg Hamiltonian, 
which are expected to determine the quantum dynamics 
in different real ID antiferromagnets, will have distinc- 
tive spectral responses in Raman scattering studies of ID 
antiferromagnets. The results summarized in Tables H] 
and [ni provide a means of obtaining useful information 
about different perturbations by comparing our results 
with future experimental results. 
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